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ABSTRACT 

The amplitude of cosmological density fluctuations, crgj has been studied and estimated by 
analysing many cosmological observations. The values of the estimates vary considerably between 
the various probes. However, different estimators probe the value of erg in different cosmological 
scales and do not take into account the nonlinear evolution of the parameter at late times. We 
show that estimates of the amplitude of cosmological density fluctuations derived from cosmic 
flows are systematically higher than those inferred at early epochs from the CMB because of 
nonlinear evolution at later times. We discuss the past and future evolution of linear and nonlinear 
perturbations, derive corrections to the value of erg and compare amplitudes after accounting for 
these differences. 

Subject headings: cosmic flows, cosmological parameters from LSS, cosmological parameters from CMBR, 
galaxy dynamics 



1. Introduction 

The as parameter is one of the most important 
and least well known parameters in cosmology. It 
is a convenient measure of the degree of the inho- 
mogeneity of the Universe. It is the rms matter 
density contrast in a sphere with a comoving ra- 
dius of 8h~^ Mpc at present, where h is the usual 
dimensionless Hubble constant in units of 100 km 
s~^ Mpc~^. The original motivation for the use 
of this parameter was the need to define the clus- 
tering amplitude: the variance in the number den- 



sity of optical galaxies in 8 h ^ Mpc spheres, cr^ 



was observed to be uni ty ( Davis fc PeeblesI 19^3t 
Strauss fc wiliid^ll99,'^^ . This parameter is sim- 
ply re lated to scale-independent linear bias (jKaiser 
19881 ). defined as the square root of the ratio of 
the number density variance to the mass variance 
b = o-gai/o-jnass = l/cTg- Later, as provided the 
standard way to normalize the theoretical power 
spectra of density fluctuations, determined from 



first principles only up to an arbitrary multiplica- 
tive constant. A theory can be tested against the 
empirical data by comparing its assumed value of 
as with a value implied by observations. A com- 
parison of values of as , estimated from various ob- 
servations can provide an important measure of 
systematic errors introduced by different estima- 
tors, based on different physical effects, e.g., bulk 
flows, cosmic microwave background (CMB) fluc- 
tuations or gravitational lensing. Such a compar- 
ison also provides an internal consistency test for 
the gravitational instability theory. 

To compare the as{z), estimated at high red- 
shift, z, with its current value, estimated from ob- 
servations at z = 0, it is necessary to take the 
growth of the fluctuations into account. Since all 
known estimates of this parameter were of the or- 
der of unity or less, until recently the linear per- 
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turbation theory was used for this purpose: 



2. Window functions 



^8(0) 



m 

D{z) 



(1) 



where D is the hnear growth factor ( Peebleslll980l) . 
As we show below, ignoring higher order terms in 
the above expression introduces a systematic error 
of order of 10 to 15%. With the latest improve- 
ments in the quality of observations, such efects 
should not be neglected and we provide a simple 
recipe how to take them into account. We focus on 
the comparison of erg estimated from the peculiar 
velocity field at an effective redshift z = 0, to erg 
implied by the CMB (effective redshift z k. 1100). 



Like many cosmological experiments, measure- 
ments of cosmic fiows are sensitive to a windowed 
integral of the matter power spectrum. In general, 
such an observable can be characterized as 



= / dkW{k)Pj,{k)T^{k) 



Relative motions of pairs of galaxies jF eldman et al, 
'200al as well as b ulk fiows (iPike fc Hudson,2005t 



(4) 



where Pp{k) gives the primordial power spectrum, 
k is the comoving wavenumber, Tf{k) gives the 
transfer function which contains the physics of the 
evolution of the particular observable from the pri- 
mordial spectrum. The window function, VF(fc), 
describes the experimental setup (sky coverage, 
depth, errors etc.). This formalism describes 



straightforward measurements of the galaxy power 

V f."^"^ ( iWatkins fc Feldman[. p^^^^.^^^ ^^^^^ ^^^^ ^ p^^^^^ 



2007t iFeldman fc Watkini l2008 t) measurements 
have been used to estimate erg and fim, the den- 
sity of the nonrelativistic matter. In particular 



from pairwise velocities we found in IFeldman et al 
|200ai 



(Tg 1.13] 



0.23 



(2) 



From CMB te mperature fluctuatio ns, the WMAP 
collaboration ( Dunklev et al. 20091 ) found 



erg =0.80 ±0.04, 



(3) 



described as the "linear theory amplitude" (see 
their Table 1). The above two estimates differ only 
slightly, at the level of 1.5-a, and it is a success of 
the model that inferences by such different meth- 
ods applied at two greatly different epochs are in 
good agreement. 

In this paper we bring the two estimates even 
closer by taking nonlinear dynamics into account. 
We find that the nonlinear correction is modestly 
significant, given the improving accuracy in cos- 
mology: the linear value of erg, hereafter denoted 
(Tl, can be smaller by ~ 10% than the nonlin- 
ear value, a systematic difference that is compara- 
ble to the current statistical uncertainties. In sec- 
tion [2] we describe the spatial window functions 
used to define a%. In section [3] we provide the 
recipe for recovering gl from a. In section 2] we 
compare our revised erg parameters to other ob- 
servational estimates of ag. We summarize our 
results in section [H 



the CMB 

spectrum for which Qi — Ci-, and the ampli- 
tude of the cosmological velocity field where Qi = 
Py{ki), the velocity power spectrum. It is cru- 
cial to note that the transfer function depends im- 
plicitly upon the other cosmological parameters 
and hence any lack of knowledge thereof will (or 
at least should) translate to increased uncertainty 
upon the amplitude. 

Bulk fiow and shear measure th e velocity- 



veloci ty power spectrurn [ or coya riance. lWatkins et al 



veloci ty DC ^ ^ 

(|2009t ): IFeldman et"aLl ' (|2009t )]. whereas pair- 
wise veloci ties measure the dens ity- velocity cross- 
spectrum (jFeldman et al.l 120031 ). Under linear 
evolution, the density contrast is proportional 
to the divergence of the peculiar velocity in real 
space, or V cx k/9 in Fourier space, so these power 
spectra differ by powers of wavenumber k from the 
density power spectrum, which can be absorbed 
into the appropriate transfer function. 

An amplitude parameter such as cg is in essence 
a spectral observable as well. We define the win- 
dow function M^jj(x, x') — W/j(x — x'), normal- 
ized so that J d^x Wr{x) — 1. For our spherical 
top hat, Wnir) = l/V where V = AttR^ /Z, when 
r = |x' — x| < -R, otherwise. Hence, the density 
contrast, spatially averaged over a sphere around 
a particular point x is simply 



d=^a;'l¥i^(x-x')(5(x') 



(5) 



The ensemble avarage of 5\ at redshift z is given 
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by 

where 



dk 
^ 



A(A:,z)2|M^(fci?)|2 , (6) 



W{kR) ^ {3/kR)ji{kR) (7) 

is the Fourier transform of the window function 
Wij(r), and ji is a spherical Bessel function of the 
first kind, while 

A = 4:TTP{k,z)k^/{2nf 



(8) 



is the dimensionless power per wavenumber oc- 
tave and P{k, z) is the power spectrum of the mass 
density fluctuations. The mean square fluctuation 
that we would measure from the actual density 
contrast depends upon the actual nonlinear power 
spectrum, but we can analogously define the linear 
variance by replacing the power spectrum in the 
above expression with the linear spectrum, con- 
strained to have the same amplitude for fc — )■ 
(i.e., on large scales where nonlinear evolution is 
negligible). 

The spatial representation of the above integral 
is given by the expression 



(9) 



function, 
(l2003h 



where ^ is the two-point correlation 
and V = AnR^/S. In iFeldman eta 
we have used the above equation to estimate 
the true present-day value of the (Ts parame- 
ter from the PSCz survey correlation function 
^Hamilton fc Teemar We have also used 

other empirical correlation functions, derived from 
different surveys, as a template and we found that 
the resulting v alue of as was not sen sitive to such 
variations; see lFeldman et al. I (l2003l) . 



3. Linear and nonlinear amplitudes 

We use two different methods to estimate the 
nonlinear corrections for erg, one based in pertur- 
bation theory, which allows us to express the cor- 
rection as a simple analytical expression, and one 
using a phenomenological mapping based on con- 
servation of pair counts and calibrated using nu- 
merical simulations, which allows us to explore the 
effects of changing many parameters individually. 



3.1. Linear perturbation theory 

Under linear evolution, the spatial and tempo- 
ral dependence of clustering separates. The den- 
sity pertuxbatim (5(x, a) = Sp/p can be described 
as (|Peebleslll980h 

Six, a) = 5'^^Hx)D{a) , (10) 

where S^^^ gives the linear density perturbation 
field as a function of comoving spatial coordinates 
X at some fiducial time and D is the growth func- 
tion, here parameterized by the scale factor a as 
a time coordinate (here and below we keep only 
the fastest-growing modes) . I n flat ACDM models 
the g rowing mode is given by ( Heatb|[l977 : Peebles! 
19931) 



Dia) 



du 



u^E^iu) ' 



where 



E{a) EE [n,na-^ + l~n,n] 



1/2 



(11) 



(12) 



In the early Universe, when the scale factor is 
small, a — > 0, equation (|lll) is well approximated 
by the expression 



D{a) = a 



(13) 



as in Einstein-de Sitter Universe. In the oppo- 
site limit, the cosmological constant becomes dy- 
namically dominant and the linear growth factor 
will saturate at a maximum value, as gravitational 
clustering is balanced by the effective force of ac- 
celerated expansion. It is easy to show that in the 
limit a — > oo, the growth factor is given by the 
expression 

2r(|)r(ii) / r!„ 

Using the above two asymptotic expressions we 
have found a new fitting formula for the growth 
factor, valid for all flat ACDM cosmological mod- 
els: 

^(«) = TTO • (1^) 



2.3 



1 + (a/^oo) 

In Figure [1] we show that equation (|15|) remains 
within two-percent level of the exact solution (ITT|) 
in both the past and the future. Note that some 
expressions for D(a) and its logarithmic deriva- 
tive, (i l og log g, frequently quoted in the lit- 
erature ( Lahav et al.lll991 : Carroll et al ]|l992l) ap- 
ply only to the past and fail for a > I. 
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3.2. Second order perturbation theory 

The perturbative solution of ttie equations of 
motion of the cosmic fluid can be written as 




Fig. 1. — The growth factor D, plotted as a func- 
tion of the scale factor a for three different values 
of ilm- The exact solution pT|) is represented by 
the solid red curve, while the dashed blue curve 
was derived from the approximate formula (|15l) . 
All models assume zero spatial curvature, JIa ~ 
1 — fljn- The scale factor at present equals unity. 
Hence the past corresponds to a/ao < 1, while the 
future to a/ao > 1 (ao = o,iz = 0) = 1). 



(16) 



where is the solution of order J and we 

assume that the linear solution is a Gaussian ran- 
dom field, so all its odd-order moments vanish. 
The mean value of the square of the above expan- 
sion is therefore given by a series of even powers 
of D(a), 



(^2) ^ ([J(l)]2)i52 + (J(2)+2^(1)J(3))^4 



(17) 



One-loop perturbative corrections to the lead- 
ing order variance gf, (r) for powB r -law power 
spectra are given by iLokas et al.l (jl996l ) and 



Scoccimarro fc Frieman ( 1996l l 



(18) 



where the factor f3 is related to the logarithmic 
slope of the two-point correlation function j{r) = 
— dln^/dlnr by the following relation: 



(3 = 1.843- I.I687. 



(19) 



The above equations and more generally, the 
weakly nonlinear gravitational instability theory 
were confirmed by N-body simulations and by 
measurements of the galaxy ske wness and bispec- 



trurn i n redshift surveys (s e e e.g.lJuszkiewicz et al 



2001 



1993: Scoccimarro et al.1 I2OOI : Feldman et al 



Verde et al.l I2002L and references therein) . 



For non-power-law spectra, 7 is a slowly varying 
function of scale. A convenient representation of 
the correlation funct ion over sc ales of interest has 
two power laws (Ha milton fc T cgmark 2002), 



(20) 



where Xj = r/rj, ri = 2.33 h Mpc, r2 — 
3.51 /i-i Mpc, 71 = 1.72, 72 = 1.28, and q = 
ct(8/i~i Mpc)/0.888. This particular choice of set 
of parameters are based on the PSCz survey, and 
this was the choice made in the original paper, de- 
terminin g from the pairwis e peculiar motions of 
galaxies (jFeldman et al.1120031 ) . We have also con- 
sidered other observational estimates of ^(r) and 
found that the resulting values of erg and 51^ were 
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unaffected. The effective slope of the correlation 
function is then 



_ dinar) 
7(r) = -- 



7ia;i 



71 



12X2 



72 



dlnr 



(21) 



independent of the amplitude q. At r = 8 h"-^ Mpc, 
the effective slope is 7 = 1.393, for which 
/? = 0.216. To map from to aj^, we invert 
equation (fT8|) : 



2/3 



(22) 



Note that for /3 — 0, the above expression gives 



= a'^ , as it should. In iFeldman et al.l 
we obtained a = 1.13. Using this value with /3 = 
0.216 in equation (|22l) . we obtain for the central 
value 

(7L = 1.02, (23) 



a modest but significant decrease. 

3.3. Phenomenological approach 

Our other method of relating linear and nonlin- 
ear variance is non-pertur bative. It uses the map - 
ping of scale proposed in [Hamilton et al.l (|l99lh , 
for which conservation of mass or pair counts re- 
lates a scale in the linear regime to a nonlinear 
"coUapsed" scale 



d{r^)[l+ar)]=R^[l + <yHR)]- (24) 



The variance is then a (nearly) universal func- 
tion 

a'{R)^g[aliri^)] . (25) 
The relation is verified and the fu nction q ide n- 



tified in numerical simulations (IHamilton et al 
I991I : IPeacock fc Doddslll996l : ISmith et al.ll2003l ). 



There is an important difference between the 
perturbative and the phenomenological calcula- 
tions. The non-perturbative mapping, given by 
equation ([25| , explicitly uses the linear expression 
for D(a) and therefore in order to derive ctl from 
a, we have to choose specific values of the cos- 
mological p arameters; for technical details, see, 
for example Peacock fc Dodd^ ( 19961 ) . This is not 
necessary for the perturbative formula ([22]) : to de- 
rive aL we only need to know cr and 7, the slope 
of the two-point correlation function at present. 
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Fig. 2. — Fully evolved variance a plotted vs. lin- 
ear variance ctl- For comparison, the isolated red 
solid line shows cr = ctl. Second order perturba- 
tion theory, based on Eq. ([T8| with /3 — 0.216 gives 
the long-dashed black line. The remaining curves 
are derived from the phenomenological mapping, 
given by Eq. The differences between them 

show how sensitive they are to different assump- 
tions about the cause of change in the amplitude 
of density perturbations. The growth induced by 
gravitational instability acting on density pertur- 
bations with the 'standard' ACDM initial condi- 
tions (rim — 0.3, = 0.7, h — 0.7) and the ap- 
propriate D[a) factor, is shown by the blue short- 
dashed line. The solid red line differs in redshift 
or amount of evolution (and hence saturates at 
0.95); the dotted black line differs in Hubble 
constant (h = 0.47-1.7) or scale {R = 5-19 Mpc) 
which changes rieff and thus 7 and so tracks a lit- 
tle differ ently from the others . We also show the 
WMAP ijPunklev et al.ll2009[ ) resuh (ctl = 0.80, 
which maps to cr = 0.85), and that from measure- 
ments of cosmological flows (cr = 1.13 correspond- 
ing to aL — 1.02, with the vertical and horizon- 
tal intervals shown in long dashed black lines the 
range in Eqs. ^ and ([271) '). 
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In Figure [2] we show nonlinear cr plotted against 
the inferred linear a (bundle of curves); for com- 
parison the isolated solid red line shows ctl. The 
perturbation theory result of equation ([^ is 
shown by the long-dashed black line. To verify 
that perturbative results still make sense even at 
the threshold of the validity of perturbation the- 
ory, when (Ti « 1, we also plot the phenomeno- 
logical mapping results, based on equation ([25]) 
for a variety of parameters. In this mapping we 
must identify why aL has changed, which may be 
from a change of fluctuation amplitude induced 
by the standard gravitational instability, evolu- 
tion epoch, or scale; each has a slightly different 
effect. In general, the perturbative curve agrees 
well with the phenomenolo gical results, based on 



the transfer function from Bardeen et al. 



and the iPeacock 
the 
the 



Doddd (|l996r ) fitting. Using 

isenstein fc Hul (119981 ) transfer function and 



Smith et al. I (l2003ir fitting make very little 



difference in the results except when the Hubble 
constant or scale changes. This occurs because a 
change in scale substantially changes the value of 7 
and so of /3. This is yet another confirmation of the 
reliability of our perturbative calculations. It also 
shows that the PSCz ^(r), assumed in the pertur- 
bative calculation, agrees well with the 'standard' 
ACDM £,{r), assumed in the phenomenological 
mapping. 

In the phenomenological ACDM relation, the 
nonlinear signal range 



maps to 



1 10+0.22 
-'^■^'-'-0.231 



^ 1 no+016 
o-L = 1.02_o^i8. 



(26) 



(27) 



Note that since a is steeper than ctl, the range in 
(Tl is somewhat narrower than the one in cr. These 
values and the resulting limits are also shown in 
Fig.H 

More generally, the existing descriptions and 
ansatzen of nonlinear clustering were created to 
describe the evolution of clustering from the past 
through to the present-day; they are not necessar- 
ily adequate representations of future clustering, 
even in the very near future. In Figure [3] we show 
the value of A^(fc, z), the mean square of the di- 
mensionless de nsity fluctuation a t various epochs 
using bo th thelPeacqck fc Doddsl ^996) mapping 
and the ISmith et al.l (|2003f ) formula; the former 



becomes time independent whereas the actual evo- 
lution will concentrate more and more matter in 
smaller haloes over time. We also see this in Fig- 
ure m which shows the value of both as and cti as 
a function of expansion factor. The latter scale 
of lh~^ Mpc is nonlinear today and therefore al- 
ready shows a difference between linear evolution 
and the different nonlinear formulae, whereas the 
canonical 8h^^ Mpc scale is just going nonlinear 
today. 




k (Mpc-') 

Fig. 3. — The linear (short-dashed black 
lines) and nonlinear values of the dimensionless 
A^(fc) as obtained from the nonlinear ansatz of 
IPeacock fc Dodds ( 1996, lo ng-dashed red lines) 
and from Smith et al. [ pOOal sohd blue lines), for 
expansion factors a = 1/4, 1/2, 1, 2, 4, and 8 
(where the present is a = 1). 

These differences arise partly from the fact 
the various fitting formulae did not attempt to 
reproduce future clustering, but moreover, from 
the well-known fact that we live at a special 
time: in ACDM, linear clustering will "saturate" 
soon. Hence, any pre diction such as that of 
Peacock fc Doddsl '|l996h based on a mapping of 



the linear spectrum will also necessarily satu- 
rate. Moreover, the various densities are cur- 
rently evolving very rapidly. For example, at 
a = 1/2, n„i was approximately what VIa is to- 
day; at a = 2, r2„i will be approximately what 
f2b is today. This also gives us some insight into 
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Fig. 4. — The value of fis (lower set of curves) 
and (Ti (upper set of curves) as a function of 
expansion f actor, for linear theory ( short-dashed 
black lines). IPeacock fc Doddsl (Il996 . long-dashed 
red lines) and from ISmith et al. (j2003i solid blue 
lines). 



the underlying (rather than just practical) lim- 
itations of these methods. The mapping pro- 
posed by iPeacock fc Dodds ( 199d ) essentially ap- 
plies equally to all scales. However, in an evolved 
ACDM universe, large-scale dynamics are dom- 
inated by the accelerated expansion and small 
scales by highly nonlinear stable clustering, which 
is captured somewhat better in the halo model. 



4. Discussion 

We will now compare different estimates of erg. 

CMB. The erg is estimated from the CMB 
flutuations, observed at z 1100 and then 
linearly extrapolated to the present era. The 
CMB itself measures a complicated (albeit lin- 
ear) functional of the power spectrum, and be- 
cause of degeneracies in the determination of 
the CMB power spectrum from the cosmologi- 
cal parameters, the value of erg from the CMB 
depends on details of the data-analysis proce- 
dure, in particular on the assumed Bayesian 
priors on the cosmological parameters consid- 
ered. The "recommended" value for a flat ACDM 



model (http://lambda.gsfc.nasa.gov) is erg = 
0.80 ±0.04. However, considering different models 
and priors on those parameters can give variations 
over 0.7 < erg < 0.9. 

Cosmic flows. Observations of peculiar ve- 
locities of galaxies provide an e stimate of erg 



at the present, nonlinear regime. [Masters et al 
(|2006() estimated crg^lj^^ = 0.52 ± 0.06 using the 
peculiar velocities fr om a sample o f clus ters in 
the SFI++ survey; IZaroubi et al.l ( 2002 ) using 
density-density and velocity- velocity comparisons 
found CTg^2°f = 0.45±0.05 both of w hich translate 
toCTg > 1. At scales of 100 Mpc. iLavaux et aL 
(|200k ) found erg = 1.72 ± 0.56 from 2MASS, al- 
though they only directly measure the velocity 
field within 30 ft." ^ Mpc, on which scales they es- 
timate a somewhat lower fluctuation amplitude. 
As mentioned above, pairwise velocity analysis 
(jFeldman et al.l 120031 ) estimates a present day 



nonlinear value of trg — 1.13lg'23- Recent bulk 
flow measurements u sing the best availab le pecu- 
liar velocity surveys (jWatkins et al. I l2009h requne 
the nonlinear erg > 1.11 (0.88) at a 95% (99%) 
confidence level. Thus, flow measurements give 
higher values of cg than the CMB, as expected 
from the ideas discussed above. We expect more 
precise measurements of erg from cosmic flows in 
the near future. The new surve ys are deeper, 
denser, have better sky coverage (jMasters et al 



I2OO6!; ISprin gob et all 2007, 2009), and we have 
improved our understanding of the distance in- 
dicat ors needed to extrac t the pecu liar veloci- 
ties ( Feldman et al. 2003t iRadburn- Smith et al 



200i iPike fc HudsonI l2005t ISarkar et al.l 12007 



Watkins fc FeldmanI f2007t iFeldman fc Watkins 
20081 ). 

Other cosmological probes. For compar- 
ison, we provide a table, where apart from the 
above estimates of erg, we also show the val- 



ues, derived from Lyo; observations (jTytler et al 
20041): cluster numb er density measurements 
2009f) : weak lensing measure - 



(Vikhlinin et al. 



ments of cosmic shear ( Benjamin et al.l 12007 ) 



the Suny aev-Zeldovich effect us ing the ACBAR 
receiver (|Reichardt et al.l l2009t) . and from the 



galaxy clustering power spectrum (ICole et a.1 
2005l : lTegmark et al'llioollEisenstein et al.lbOOsI^ 



All of these measurements fall in the range 
0.7 < erg < 1.2 with the most likely values in 
the interval 0.8 < erg < 0.95. 



7 



Table 1: as from various estimators 



Method 


Parameter 


value 


CMB^ 


o-L 


0.80 ±0.04 


LY-a"! 


c^L 


0.85 ±0.02 


Cosmic Shear ^ 




0.84 ±0.05 


Chisters ^ 




0.75 ±0.01 


SZ (ACBAR) ^ 


0-8 


n Q4+0 03 
^■^^-0.04 


Galaxies ^ 


gal 


0.92 ±0.06 


Flows 






Pairwise velocities J_ 


0-8 


-1 iq + 0.22 
^•^"J-0.23 


Bulk flow 


0-8 


> 1.11 (0.88) at 95% (99%) CL 



(|2009[ ) 



Dunkle v et al] |200£ ); ^Tvtler et al.l (|2004|): teen iamin e t all (l2007l): Wikhlinin et all (l2009l) 
Reicha rdt et all I2OOS h ICole et al.l (|2005l ) nibgmark et al.. (2004 ) : lEisenstein et all (|2005[ ): iFeldman et al 
(j2003h : ^Watkins et al 



5. Conclusions 

We have presented a formalism to calculate as , 
the amplitude of cosmological density fluctuations 
on scales of 8 Mpc, that incorporates the non- 
linear evolution of the parameter. Estimates of 
as depend directly on the epoch and scale of the 
surveys used. When using deep, high-redshift sur- 
veys (CMB) that estimate as in the linear regime, 
the results suggest systematically lower values of 
as- When analyzing shallow, local data (pecu- 
liar velocities) which are affected by nonlinearities, 
we get higher erg. The results from other cosmo- 
logical probes, shown in Table 1, show a similar 
trend. When results from deep surveys are being 
corrected for this effect, most estimates from vari- 
ous independent surveys on all scales agree better 
with each other: the systematic differences are re- 
duced. 

Quantitatively, our main result here is the re- 
duction in the value of as, derived from the ob- 
served mean pairwise velocity of galaxies. The 
original estimate, as we have discussed earlier, is 

as = 1.13±EJ:1 . (28) 

After the correction based on the second order per- 
turbation theory, this becomes 

aL - imt^oil (29) 

Somewhat more cumbersome non-perturbative 
methods give identical results, bringing late-time 



estimates of crs close to high-redsift measurements, 
which appear in Table 1. 
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